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An effect of geometrical pliase sliift is predicted for a light beam propagating in the field 
of a gravitational wave. Gravitational radiation detection experiments are proposed using 
this new effect, the corresponding estimates being given. 
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^ ; 1 Introduction 



O 

"q I M. Berry (1984) has shown that a quantum system whose Hamiltonian H{() depends on some 

O^' parameters C"" and which evolves in time in such a way that during the evolution the state of the 

5—1 . system traces out a closed curve C in the space of these parameters, the wave function can get 

an additional geometrical phase 0{C). This geometric phase depends on motion of the system 
in the space of parameters. 

^ . Recently an effect of rotation of the polarization vector was predicted (Chiao and Wu 1986) 

^ \ for a linearly polarized laser beam travelling along a single helically wound optical fiber, the 

results having been immediately experimentally confirmed and their connection with Berry's 
phase (Tomita and Chiao 1986) shown. Later it was found (Cai et al 1990) that this Berry's 
phase has in fact a classical origin, and it arises from the intrinsic topological structure of 
Maxwell's theory if the Minkowski space-time is considered as a background. If k is the wave 
vector of an electromagnetic wave and e(k) its complex polarization vector, then the condition 
k^ = u;^(k) = const determines a sphere S^. The angle ^(k) of rotation of e(k), when it is 
parallelly transported along a null geodesic with the tangent four-vector A;, is defined by 

56 = ieL(k) • deij(k), 

where 

gr = 2~^/^(ei + ie2), ei = 2"^/^(ei - 262), x eR = ik/u 

and e = e~'^e^. The integral angle occurs to be = j 56 where integration is performed over 
the path of k on S^; thus 6 has a geometric origin. When the curve is a closed path C on S*^, the 
angle 6{C) is given by 6{C) = — r2(C), where Vl{C) is the solid angle of the loop C with respect 
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to the center of the sphere. The expression above is essentiaUy a flat space-time expression (see 
(Cai et al 1990)), however apphcable also to general relativity (cf. (Bildhauer 1990)). 

In the second example of the Berry's phase in optical experiments (Pancharatnam 1956), the 
state space is the Poincare sphere which describes all possible polarization states of light. For 
this case the direction of the light propagation is fixed and a cycle of changes in polarization 
states corresponds to a closed curve on the Poincare sphere (see e.g. (Bhandari and Manuel 
1988, Simon et al 1988)). The phase 9{C) is known as Pancharatnam phase (Pancharatnam 
1956) and given by e{C) = -{l/2)n{C). 

Here we predict a similar effect of geometrical phase shift for light beams propagating in 
the field of a gravitational plane wave or pulse of gravitational radiation: 1) cyclically from 
and (after reflection) to an observer thus being closely related to the Pancharatnam phase; 2) 
along a circular flbcr of radius Rq. Wc show that for a light beam orthogonal to the direction of 
propagation of a gravitational wave, this phase grows proportionally to L/X where L is distance 
between the observer and reflecting system and A, characteristic wavelength of the gravitational 
wave packet. In the second case if A = 7ri?o, the resonance occurs and the phase shift grows 
proportionally to m, the number of revolutions of light. For a pulse of gravitational radiation, 
the relative phase shift is proportional to the characteristic amplitude of the pulse in either of 
these two cases (for preliminary results see (Mitskievich and Nesterov 1995)). 

In this paper we use the space-time signature (-1-1, —1, —1, —1); Greek indices run from to 
3 and Latin, from 1 to 3; it is essential to remember these notations when the integration by 
parts is performed (see (4) and, in the Appendix, (9)). 

2 General results 

A concise description of this phenomenon can be done using the Ncwman-Pcnrosc formalism 
which is applicable to propagation of light in arbitrary media (we are interested here in the 
cases of (1) a vacuum and (2) an optical flbre, the both in a gravitational field). The real null 
Newman— Penrose vectors arc I = (k^, k) (tangent to the light world line) and n {I ■ n = 1), the 
complex ones, m and fn (m ■m = —1). Wc shall consider m = e~^^eji to be the circular (right) 
polarization vector. When light propagates along optical fibers (not necessarily geodesically), 
the Newman-Penrose description is quite essential for differentiation along I, D = Vi being a 
Newman-Penrose operator. The corresponding equations reduce to 



Coefficients in the right-hand side of these equations are in general complex functions (a bar 
denoting complex conjugation), e — e determining torsion and n curvature of the space-like 
trajectory of light (an optical fiber); cf. (Penrose and Rindler 1986, p. 169 ff). The change of 
polarization angle generally reads 



Dl = (e + t)l — Rm — Kfh, 
Dn = — (e -|- e)n + vrm + 7fm, 
Dm = ttI — Kfi + {e — e)m. 
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In a vacuum, a light beam propagates geodesically, thus Dl = 0, and equation of the polar- 
ization vector transport reads Dm = (cf. [10]); these two equations yield also Dn = 0. For a 
planar optical fiber, Dm = itl — nn. Hence in the both cases fh ■ Dm=0, thus 



6* = iy eL • DeRdr], (1) 

the integration being performed along the light world line T canonically parametrized by rj. For 
the left polarization one has to exchange subscripts L and R in (1), or, equivalently, to change 
the sign in the right-hand side of this equation. 

Note that under the parallel transport along the space-like geodesies orthogonal to the ob- 
server's world line 7, the vector of right (left) polarization (e^) does not change, but it 
changes under the transport along T (the null line of light whose characteristics arc measured 
in course of the proposed experiment). This very fact makes the existence of nontrivial phase 9 
essential. 

If we intend to consider an experiment of detection of gravitational radiation, it is natural 

to connect the Newman-Penrose frame with the Fermi coordinates X^^ = 6qS + S^^^u (see e.g. 
(Manasse and Misner 1963, Misner etal 1973)). Here s is proper time along the observer's 
geodesic world line 7, u being proper length parametrizing the (space-like) geodesic orthogonal 
to it, with a unit tangent vector ^ on 7. In fact, ^ describes the direction in which such a space- 
like geodesic goes, and it possesses only spatial components different from zero (the temporal 
coordinate = T is directed along 7). In Fermi coordinates, components of the corresponding 
orthonormal tetrad e(y^^ parallelly transported along the spacelike geodesic, are represented as 
expansions 

H^f = ^^-\ ^'"ijo X'X^ + ■■■, 

and the connection coefficients (Christoffel symbols) take the form 



^to =R^viO X'^ -\ , 

rf, = -^^%> x' + ---. 



Quantities of the type of Q are taken on 7. 

The radius of convergence of series is determined by conditions (Manasse and Misner 1963) 



uq <C min < 





1 I Ra-ySP 

RajSp I ' I Ra-yS/Sji 



(2) 



, 

The first condition uq <^ \ Ra-ySfS \ ' determines the size of V where the curvature has not yet 
caused spatial geodesies to cross each other. The second condition defines the domain where the 
curvature does not change essentially. For instance, for gravitational waves with wavelength 77 
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the Riemann tensor is ~ ^exp(iA:^a;'^)/A^, where A is the dimensionless amphtude. So equation 
(2) yields 

uo < mm{X/VA, A}. 

Generally it is assumed A < 10~^^. This means that the size of V is restricted by uq <^ A. So the 
application of Fermi coordinates to the modern experiments may be very restrictive since A is 
often supposed being in the order of 300 km. Thus for enlarging the range of validity by a factor 
1 / y/A (which is about 10^ in our example) it is necesary to take into account all derivatives of 
the Riemann tensor. 

Covariant derivatives of the Fermi basis can be described as (see the Appendix) 

^d.e^.f = - r R^uXp^dr + - Tdr [ ddta^^R^' ^i^e dr' + 0{R^) (3) 

JQ U Jo Jo 

where integration is performed along a space-like geodesic orthogonal to the Fermi observer 
world line 7. This integral formula obviously includes all derivatives of the Riemann tensor. We 
shall apply equation (3) for calculation of the phase shift (1). 
Using (3) we find from (1) 

9 = -i [ CL^en'^l^dX f R^^xpedr + i [ -eL^CR^dri ^ dr [ R^^i/^^dr' + 0{R^). 
Jr Jo Jr Jo Jo 

Integrating by parts we obtain 

9 = -i [ BL^CR'^fdrj [ Ri,,opedT + i [ -eL^eR^l'dr] [ R^^ip^rdT + 0{R^). (4) 
Jr Jo Jr Jo 

We consider now a plane weak gravitational wave whose metric tensor is usually written in 
synchronous coordinates, 

ds"^ = 'q^ydx^dx'' + hahdx"-dx^, 
where a and h run from 1 to 2 while 

Kb = habit - z), h22 = - fill = h+, hi2 = hx- 
Using the definition of the Riemann tensor 

we find that in the linearized theory non-zero components of Riemann tensor are 

RsabS = RoabO = —RsabO = ^Kb: Rti22v = —Rfxllv, R/xUiy = RfjZlv (5) 

where dot being derivative with respect to t. Now the equation (4) is readily applicable, and two 
typical cases emerge: (A) parallelly (antiparallelly) propagating gravitational and light waves, 
(B) mutually orthogonal waves. Below we study the both cases. 
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2.1 Parallelly (antiparallelly) propagating gravitational and light waves 

We assume that the both waves propagate along z axis. One can write / = /^(l, 0, 0, ±1), where 
Z° = dT /drj; the upper sign corresponds to positive direction propagation of hght and the lower 
sign, to negative one. Let us take 

6^/^ = -^(0,l,i,0), ei'^ = -^(0,l,-i,0). 
Then equation (4) reduces to 

e = - dr R^20i,edr'± — - / i?i236eVdr' + 0(i?2). (6) 

Jo Jo Jo u{t) Jo 

Applying (5) we find 9 = that means absence of the phase shift. 

2.2 Orthogonally propagating gravitational and light waves 

Let the gravitational wave propagate along z axis and the light beam, in the (x, y) plane at an 
angle (f) to x axis. We suppose 

^ = (0, cos 0, sine/), 0), / = /°(l,±cos(^,±sin0,O), f = dT/dri, (7) 

1 1 
eij^ = — p(0, ^sin^, ±cos^,z), e^'^ = — ^(0, ^ sin ^, it cos 0, — i), (8) 
V 2 V 2 

where the upper sign corresponds to propagation of light from the observer, and the lower one, 
to the observer. Then the phase shift (4) is given by 

j-T i-u{t) 

e = ± dr {Rsuoicos'^ (t> - sin^ (f>) +2R3220 sin 4) cos (f>) dr' + 0{R^) (9) 
Jo Jo 

where we used (5). It is convenient to rewrite this equation as 

e = ^- dr (/ixCOs20 + /t+sin20)dr + O(i?2). (10) 
2 Jo Jo 

Taking into account that X'^ = a;'* + 0{h) and h = h{t — z) {h being hx or hx), one can write 
(10) as 

6 = ^]- f (hxiT)cos2(f) + h+{T)sm2(j))u{T)dT + 0{R^). (11) 

Let us define 6 = a_|-sin2(/) + ax cos 2(/), where a+ and correspond to the two independent 
polarization modes of gravitational wave (see e.g. Misner et ctl 1973). Then (for ot^ or ctx) 

a = Tl£h{r)u{T)dT + OiR^), (12) 



h being h+ or hx- Integration by parts with du = ±dT yields 

T 




a = \ {h{r) ¥ /i(t)u(t)) 
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where, as above, the upper sign corresponds to propagation of (right polarized) hght from the 
observer, and the lower one, to the observer. For the left polarized light the overall sign of the 
right-hand side of the both last equations should be changed. 

Let us consider an experiment when a circularly polarized electromagnetic wave is travelling 
between the observer and a reflecting system in the {x, y) plane, the gravitational wave propa- 
gating in the positive direction of z axis. If (say, right) polarization does not change in course 
of reflection, the phase shift is 



Aai = ^[h{2T)-h{0)-2Th{T)], 



and if polarization changes to the left one, 

1 



Aa2 = 7:[2h{T) - h{2T) - h{0)] 



(13) 



(14) 



(the effects for an initially left polarized beam have inverse signs). 
For multiple reflections we have 



1 ^ 
Aai(Ar) = -[h{{N + l)r) + h{NT) - h{T) - h{0) - 2r ^ h{mT)] 



m=l 



Aa2{N) = -HNT) - h{{N + 1)T) + h{T) - h{0)], 



(15) 
(16) 



where N is the number of reflections. 

For a plane monochromatic gravitational wave, h{t — z) = Acos{Lo(t — z) + S) where A is a 
dimensionless amplitude. In this case (15), (16) are rewritten as 



Aai(iV) = ^(A-sinA) 



sin(f A)sin(^A + (5) 
sin(A/2) ' 
N + 1 

Aa2iN) = 2Asin(A/2)sin(ArA/2)cos(^^A + (5). 



(17) 
(18) 



Here A = 2nL/X, L being distance between the observer and reflecting system, and A the 
gravitational wavelength. 

The average relative phase shift between the right and left polarized light, is 



y^< AOi^ > = A 



(A - sin A) 



sm 



sin(A/2) 



2A 



A NA 
sm — sm — — 
2 2 



(19) 
(20) 



A = y A'^ + A^. being the dimensionless amplitude of an unpolarizcd gravitational wave, and 
averaging is performed with respect to all polarizations and phases 5+ and 6x ■ 

We shall consider here no phase change in course of reflection since this is the only way to 
obtain a sensible intergal effect. If A S> 1, then, for instance for A?^ = 1, we obtain 



V< Ae^ > = AA. 



(21) 
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We see that effectively the dimensionless amphtude A grows by a factor A. When N ^ 1, we 
find that maximum of phase shift 



V< A6'2 > = 2TTmNA 
occurs for A = 27rm where m is integer. If A ^ 1, we obtain 



V< y = 2A 
and relative shift takes its maximum value, 



sm 



V< > ^ 2A, 



(22) 



(23) 



(24) 



when NA = ir. 

Another possible experiment involves a scheme similar to that of (Braginski and Menskii 
1971) but with measurement of the geometric phase shifts for right and left polarizations (and 
not the frequency shift as in [15]) of light travelling along a circular fiber of radius Rq. Let us 
consider the case when a plane gravitational wave is propagating in the positive direction of z 
axis, the fiber lies in {x,y) plane and the light propagates in counter-clockwise direction. We 
suppose 



^ = (O,cos(?!),sin(?!),0), Z = /°(1, - sin<;/i, cos ^, 0), f = dT/dri, 



en'' 



1 

72 



(0, cos0,sin^, i), ei,'* 



V2 



(O,cos0,sin0, — i). 



Then the phase shift is given by 
rT 



e 



Jo 



R3awe{r)e{r)dT + 0{R' 



Ro 



/o 2 jQ 

where we used equation (5). Defining 9 = + Ox we obtain 



(25) 
(26) 



(27) 



a^ = ^ I h+{T) cos (2a;or)dr + 0{R^); 
2 Jo 

ax=-^^ h^{T)sm{2L00T)dT + O{R^), 

where uq = 27r/To, and Tq is the period of revolution of light. Integration by parts with 
/i(r) = Acos{uT + S) where, as above, A is dimensionless amplitude, yields 



A, 



■sin((a;o - uj/2)T) cos((wo - oj/2)T - 6+) 
1 - 2{uJo/uj) 
sin((wo + <^/2)T) cos((ljo + w/2)T - 6+ 
1 + 2{loo/uj) 

( ^ \ r sin((a;o + lo/2)T) sin((a;o + uj/2)T + ) 
V2a;oy [ 1 + 2(ljo/cj) 

_ sin((a;o - uj/2)T) sin((a;o - uj/2)T - Sy, 
1 - 2(a;o/w) 



(28) 



(29) 
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For left polarized light, the overall sign of the right-hand side of the both last equations should 
be changed. The average relative phase shift between right and left polarized light, is 



A uj 



2^2 



siv? 



/2)r 



+ 2 



sin(a;o - a;/2)rsin(u;o + a;/2)rcos(2a;T) - uj/2)T 



.1/2 



+ 



(1 - 2(a;o/a;))2 J 



(30) 



When the gravitational wavelength is equal to 7ri2o, a resonance occurs leading to V < A^^ > 
\/2 irmA where m is the nTimbcr of revolutions of light. It is clear that such a detector is effective 
for high frequency gravitational radiation, the corresponding Rq being around 100 km for 10'^ 
Hz. A considerably smaller size of detector could be achieved for a toroidol winding of the fiber, 
this case being currently under consideration. However measurements of the phase shift in a 
fiber (and in any media other than a good vacuum) could be virtually impossible due to the 
random fluctuations, so that we consider here these cases to the end of completeness only. 



3 Discussions and conclusions 

We would like to discuss the possibility of detecting gravitational radiation based on the proposed 
new effect. Let us consider the case A ^ 1. From (21) we obtain V< A6'2 > 2 • m~^ALv, 
where v is characteristic frequency of the gravitational radiation in Hz and L is the distance 
in km. It is clear that such a detector is effective for high frequency gravitational radiation 
{y ~ 10^ Hz). If, for instance, the reflecting system is placed on the surface of moon, we 
have V < A02 > ^ ^Av. For the 5-million- kilometer- long Laser Interferometer Space Antenna 
(LISA), which would fly in heliocentric orbit (see Thorn 1995a, b), our estimation of phase shift 
is V<Ae^ ~ 10^ Av. 

Let us compare the experiment proposed here, with the standard experiments involving 
resonant antennae directed to the Virgo cluster and tuned to some 3000 Hz. In this case 
A ~ lO"^'' (see the corresponding data in (Douglas and Braginsky 1979, Thorne 1995a,b)). So 
the geometrical phase detector with the base Earth - Moon treats this radiation as if it had 
an effective magnitude of some 10~^^ and for LISA as if it had an effective magnitude of some 
10-1^ 

For experiments using a laboratory size apparatus, A <C 1, we see that V < AO"^ > ^ 
2A |sin(A/"A/2)|, and the relative shift takes its maximum value when NA = tt. For LIGO/VIRGO 
(L ~ 3 • 10^ cm) interferometers we find that the correspondent frequency of the gravitational 
wave is ~ 3 • W^/N . It is known that for L « 10^ cm one could expect about 300 reflections 
(Douglas and Braginsky 1979) which is sufficient for detection of continuous waves with 100 
Hz < < 10^ Hz. The advanced LIGO interferometers are expected to have their optimal 
sensitivity at z> ~ 100 Hz, and rather good sensitivity all the way from v ^ 10 Hz to u ^ 500 
Hz (Thorne 1995a,b). 

If we are interested in detection of a burst of gravitational radiation with L <C A (let its 
maximum be at the moment of time NT/2 and characteristic duration r = NT), then (15) 
yields V< M"^ > ~ A for N ^ 1 or N = 1, where A is characteristic dimensionless amplitude 
of the pulse. Realistic estimates for L and (LIGO/VIRGO) are L 3- 10^ cm, W<N < 300, 
while the characteristic frequency is 100 Hz < v < 10^ Hz. 
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Measuring of this effect consists of a comparison of interference patterns for the both circulary 
polarizations of a light beam. Similar measurements of a phase shift between opposite senses of 
circularly polarized light were performed using a nonplanar Mach-Zehnder interferometer (Chiao 
et al 1988); a very clear exposition of theoretical and experimental details concerning observation 
of phase shifts due to geometrical and topological effects, see in (Chiao 1990), naturally, without 
dealing with any gravitational effects. 

It is worth making a comment on the order of smallness of the dimcnsionlcss amplitude of 
the gravitational wave A in connection with the limits of short gravitational waves. This order 
is most invariantly attributed to the space-time curvature which is calculated using the metric 
tensor, see (5): Riem ~ A/X'^. Thus when A — ^ 0, A should also tend to zero if the curvature 
keeps its order of magnitude unchanged, and in the limit of short gravitational wavelengths 
(geometric optics in the gravitational sense) the apparent divergence of the predicted effect at 
small A's (see equations (17), (21)), is merely spurious: 

The effect we predict in this paper, makes it in principle possible to detect gravitational waves 
using not an interferometer as a whole, but only one of its arms, since there is a fundamental 
difference in propagation of the left and right polarizations along one and the same null line of the 
light. One has merely to separate the light of these different polarizations after it has returned 
from its travel, then to transform the (circular) polarization of one of the resulting beams to the 
opposite one, and finally to observe the interference fringes after mixing the beams. 

We think that this effects reflects an interaction between photon's spin and the space-time 
curvature, which is closely related to the well-known Papapetrou-Mathisson effect. 
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Appendix 

Here we shall obtain the formula (equation (3) in the text) 

Va,e(,)'^ = - r Rt^,^^edT + - r dr f dr' SlR^^^i^e + 0{R^). 
Jo u Jq Jo 

Let us start with the Taylor expansion for tetrad in a world tube surrounding the world line 
7 of inertial Fermi observer: 

d e/.A 1 d ei,,\ „ \ or ei,,\ „ 

u being the canonical parameter along spacelike geodesies orthogonal to 7. Using in Fermi 
coordinates the equation of parallel propagation 

^ + r^.e(.)T = 0, (32) 

we rewrite equation (1) as 



9 



The expansion of the connection coefficients is given by 



r;;^ =Kx,i + ^ x'x' + .... (34) 



Applying equations (3), (4) we obtain the following series for the spatial covariant derivatives 
of tetrad 

^0 jtl° ° k I 

+ l^{r>^u^,k,l,m - r^(i,fe,;,^))X^X^X™ + 0{R'). (35) 

Now using the relations 

2 ^'^ vkiX = {T^ui,k - ^'^u(i,k))X , 

2 , , , , 

DW vk vl /"p/i T^ti \ vk Yi 

2 ^ vki,l^ — V- vi,k,l — v{i,k,l))-^ 

30 k I m ^ ^ k I m 

^ R'^ uki,l,mX X X™ = {T'^ui,k,l,m - u{i,k,l,m))X X X"^ , etC 

one can write the expansion (5) as 

I 2 30 
'^di^{v)'^ = 2 vkiX + ^ -^'^ vki,lX X + — R'^ uki,lX X X"^ H 

+ (;7Tl)! ^''-h^M,■■■,lnX''X''■■■X'-+■■■ + 0{R'). (36) 

It is convenient to present this series in the form 





^^^'M = \ di? (n + 2)(n!) + ^ 

n=0 ^ / \ / 

Straightforward calculation yields the following integral representation of equation (7): 

1 

Va,e(,)^ = — / R'^uiki'^TdT + 0{R^). 

Integrating by part we find 

^d,e^uf = - r R^uik^^dr + - TdT r dr'R'^.ike + OiR"). 
Jo u Jo Jo 



(37) 



(38) 



(39) 



Now let us calculate the temporal covariant derivative V^ge^,^)'^. From equations (2)-(4) we 
easily obtain 

fc 1 ° ° k I 

"^do^iiy)^ =^^uO,k X + ■^{T^uO,k,l - ^^uk,0,l))X X 

+^(r^o,fe,;,m - r^fe,o,;,m)^'^'x™ + • • • + +oiR^), (4o) 
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where we have taken into account that r^,^o= 0. Obviously, 

+^ R"" .ho,h,- ,inX^'X^^ . . . X'" + • • • + C»(i?2), (41) 

or 

n=0 ^ ' 

The correspondent integral representation reads 

^doH-f = - r R^-ok^^dr + 0(i?2). (43) 

Finally one can combine equations (9), (13) and write 

^d,e^.f = - r R^'uxpedr + - Tdr f dT'6\R'',ipe + OiR"). (44) 
Jo u Jo Jo 

Noting that the transformation from arbitrary coordinate system to the Fermi one takes the 
form = A^^X'^ + 0{T) (or = ^^u + 0{T) ) one can write Eq(14) in an arbitrary coordinate 
system as 

^d,e^.f = - r R\sp{^-'raA.\A.\edr 
Jo 

+- r dT r dT'R%ipiA-'r^A\A\e + 0{R^). (45) 
u Jo Jo 



11 



References 

Berry M V 1984 Proc. Roy. Soc. A 392 45 

Bhandari R and Manuel J 1988 Phys. Rev. Lett. 60 1210 

Bildhauer S 1990 Class. Quantum Gravity 1 2367 

Braginsky V B and Menskii M B 1971 Sov. Phys.— JETP Lett. 13, 417 

Cai Y Q, Papini G and Wood W R 1990 On Berry's Phase for Photons and Topology in 
Maxwell's Theory (University of Regina Preprint) 

Chiao R Y and Wu J S 1986 Phys. Rev. Lett. 51 933 

Chiao R Y, Antaramian A, Ganga K M, Jiao H, Wilkinson S R and Nathel H 1988 Phys. 
Rev. Lett. 60, 1214 

Douglass D H and Braginsky V B 1979 Gravitational-radiation experiments, in General 
Relativity. An Einstein Gentenary Survey,e6ited by S W Hawking and W Israel (Cambridge 
University Press, Cambridge). 

Manasse F K and Misner C W 1963 J. Math. Phys. 4 735 

Misner C W, Thorne K S and Wheeler J A 1973 Gravitation (W.H. Freeman, San Francisco). 

Mitskievich N V and Nesterov A I 1995 Gen. Relat. Grav. 27 361 

Pancharatnam S 1956 Proc. Ind. Acad. Sci. Ser. A 44 247 

Simon R D, Kimble H J and Sudarshan E C 1988 Phys. Rev. Lett. 61 19 

Penrose R and Rindler W 1986 Spinors and Space-time, Vol. 2 (Cambridge University 
Press, Cambridge) 




Tomita A and Chiao R Y 1986 Phys. Rev. Lett. 51 937 



12 



